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duals of QCD from experimental hadron properties. This framework allows us to derive 
a gravity dual which reproduces the empirically found linear square-mass trajectories of 
universal slope for radially and orbitally excited hadrons. Conformal symmetry breaking 
in the bulk is exclusively due to infrared deformations of the anti-de Sitter metric and 
governed by one free mass scale proportional to Aqcd- The resulting background geometry 
exhibits dual signatures of confinement and provides the first examples of holographically 
generated linear trajectories in the baryon sector. The predictions for the light hadron 
spectrum include new relations between trajectory slopes and ground state masses and are 
in good overall agreement with experiment. 
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1. Introduction 

Since the inception of Quantum Chromodynamics (QCD), progress in understanding its low 
energy realm was hampered by the scarcity of adequate techniques for handling strongly 
coupled Yang- Mills theories analytically. The discovery of the AdS/CFT correspondence 
m, ^ has given promise for this situation to improve in a qualitative way. Indeed, the 
ensuing dualities explicitly relate gauge theories at strong coupling to physically equivalent 
string theories in ten-dimensional spacetimes which become tractable at least in the weak 
(string) coupling and curvature limits. These dualities are manifestations of the holographic 
principle Q and have triggered an entirely new way of thinking about nonperturbative 
QCD. 

The currently best understood dualities deal with supersymmetric and conformal gauge 
theories, however, and so far the existence of an exact QCD dual and its explicit form have 
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not been established ab initio. The pioneering apphcations to QCD and hadron physics^ 
therefore rehed on a minimal infrared (IR) deformation of the anti-de Sitter (AdS) metric to 
model confinement |^. By restricting the fifth dimension to a compact interval, this "hard 
IR wall" geometry softly breaks conformal symmetry in a way consistent with high-energy 
QCD phenomenology (including counting rules for exclusive scattering amplitudes etc.) 
Q . The extension of this minimal approach into a search program for the holographic dual 
of QCD, guided by experimental information from the gauge theory side, is often referred 
to as AdS/QCD. 

At present, this program is being pursued along two complementary lines. The first 
one is mostly bottom-up: it assumes a five-dimensional, local effective field theory in IR- 
deformed AdSs spacetime (and potentially in additional background fields of stringy origin) 
to describe the gravity dual of QCD and attempts to constrain its form and parameters 
by experimental information. For various works in this direction see ^, ^, ^, 10, 11, 12| 



and references therein. The main virtue of this approach lies in gathering and organiz- 
ing information on holographic QCD by using the wealth of accurate experimental data 
available on the gauge theory side. The second type of approach is more directly guided 
by the underlying, ten-dimensional brane anatomy of the gravity dual and attempts to 
maintain closer ties to it^. Up to now, however, this typically comes at the price of less 



direct relations to QCD. For recent work along these lines see for example |14, 15, 16 1 and 
references therein. 

Over the last years AdS/QCD has met with considerable success in describing hadron 



properties, the heavy quark potential jl^, [l^, |lyl, vacuum condensates ||2^, QCD scatter- 



ing amplitudes at high energy |^ |T5| etc. Many of these results were obtained on the basis 
of the minimal hard wall implementation of IR effects. Although very useful in several 
respects, the hard wall is also an oversimplification and has revealed shortcomings when 
more complex and quantitative QCD properties are considered. The perhaps most impor- 
tant limitation discovered so far is that it predicts quadratic instead of linear square-mass 
trajectories both as a function of spin and radial excitation quantum numbers |l^, in 
contrast to experimental data and semiclassical string model arguments (for highly excited 
states) which relate linear trajectories to linear quark confinement p^ . 

Different ways to overcome this problem in the meson sector have recently been pro- 



posed in Refs. |11, |T^, 23|. As a result, linear Regge trajectories cx J for spin- J 
excitations or the analogous radial excitation trajectories oc of mesons could be re- 
produced by different holographic models. Although the baryon sector 0, 24, 25| exhibits 
similarly pronounced empirical trajectories of the above type [^], however, a dual descrip- 
tion for them has not yet been found. Our primary goals in the present paper will therefore 
be to understand how linear baryon trajectories can arise in the AdS/QCD framework, and 
to establish dual descriptions for additional spectral signatures of linear confinement. 



^Applications to remoter "cousins" of QCD, and in particular to their glueball spectrum, have a longer 
history. See for example Refs. Q. 

^At present there is not enough calculational control over string theory in curved spacetimes to allow 
for a strict top-down approach. Nevertheless, on the maximally symmetric AdSs x a complete solution 
of the world-sheet string theory may ultimately become feasible even in the high-curvature regime 
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To this end, we will focus on another striking and systematic feature in the light hadron 
spectrum, namely the combined linear square-mass trajectories 



M^ = M^ + W{N + L) (1.1) 

of Regge type on which radial (A^) and orbital angular momentum (L) excitations join. 
The trajectory structure (^]^) is experimentally well established for both the light meson 
and baryon resonances. Two fits to the meson data yield mutually consistent mean slopes 
W = 1.25 ± 0.15 GeV2 ^ and 1.14 ± 0.013 GeV^ [||], respectively. The fit to the light 
baryon resonances (i.e. to those consisting of up, down and strange quarks) results in the 
somewhat smaller but still compatible slope W = 1.081 it 0.035 GeV^ p6| . Hence the 
value W ^ 1.1 GeV'^ is approximately universal^ for all trajectories |2S]. The ground state 
masses Mq, on the other hand, are channel dependent. 

Our work will rely on the dual representation of hadronic states with higher intrinsic 
angular momenta by metric fiuctuations [Q, |l^, |2l| which suggests itself for our purposes 
because it provides direct access to orbital excitations. The main strategy will be to de- 
vise and apply a method for deriving IR deformed gravity duals which incorporate the 
experimental information contained in the hadron trajectories ( |1.1D . After a brief sum- 
mary of pertinent facts and results from AdS/QCD and the hard- wall metric in Sec. |2|, 
we set out heuristically in Sec. by determining a minimal modification of the AdS mode 



dynamics which generates the linear trajectories (1.1) in both meson and baryon spec- 
tra. Subsequently, we explicitly construct the IR deformations of AdSs which encode the 
same dynamics, by deriving and solving differential equations for the conformal symmetry 
breaking part of the warp factor in Sec. ^. Several new features of the resulting gravity 
background are discussed in Sec. |. In Sec. I we determine the value of the conformal 
breaking scale and compare the results of our holographic description to experimental data. 
In Sec. 0, finally, we conclude with a summary of our findings and mention a few avenues 
for future improvements and applications. 



2. AdS/CFT correspondence and hadron spectrum 

The gauge/string duality |[l|, |2| maps type IIB string theories in curved, ten-dimensional 
spacetimes into gauges theories which live on the (flat) 3+1 dimensional boundaries. For 
an UV-conformal gauge theory like QCD, the dual string spacetime is the product of 
a five-dimensional non-compact part which asymptotically (i.e. close to the boundary) 
approaches the anti-de Sitter space AdSs (R) of curvature radius R, and a five-dimensional 
compact Einstein space X5 (where X5 = (R) for the maximally supersymmetric gauge 
theory) with the same intrinsic size scale. The line element therefore takes the form 

= e2^(")^ {rj^^dxf'dx'' - dz"^) + R^ds% (2.1) 

^This may be an indication for the conformal symmetry breaking scale oc Aqcd to be approximately 
hadron independent in the light flavor sector, as noted for orbital excitations in Ref. |[7|. 
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(in conformal Poincare coordinates) where 77^^ is the four-dimensional Minkowski metric. 
Since A ^ breaks conformal invariance explicitly, one has to require A{z) ^ as z — > 
in order to reproduce the conformal behavior of asymptotically free gauge theories at high 
energies. The string modes^ (p^ {x, z) = e~*^'^/j (z) dual to physical states of the gauge 
theory are particular solutions of the wave equations^ in the geometry (2.1) and fluctuations 



around it, and potentially in additional background fields of stringy origin |2[. 

Casting these wave equations into the equivalent form of Sturm-Liouville type eigen- 
value problems, one finds 

[-d! + Vm {z)] VP (z) = Ml^ {z) (2.2) 

for the normalizable string modes (-z) = g (z) Jm {z) dual to spin-0 (M = S) and spin-1 
(M = V) mesons as well as from the iterated Dirac and Rarita-Schwinger equations 

[-dl + Vb,± (z)] V'i {z) = Mli^^ (z) (2.3) 

for the string modes V'± {z) = h (z) /b,± {z) dual to spin-1/2 and 3/2 baryons (where it 
denote the two chiralities of the fermions with ij^il^^ = 0. The potentials Vm,b 

contain all relevant information on the string mode masses and the background metric ( |2.1] ) 
which also determines the functions 5 (z) , h (z) introduced above. The eigenvalues ^ 
constitute the mass spectrum of the four-dimensional gauge theory on the AdS boundary. 
The boundary conditions for the eigensolutions if and V'i are supplied by specifying the 
corresponding gauge theory operator according to the AdS/CFT correspondence Q (cf. 
Eq. (p.4|)) and by the requirement of normalizability (and minimal string action in case 
of ambiguities) of the eigenmodes. In some cases a further, less well determined boundary 
condition is imposed in the infrared, at z = Zm, in order to break conformal symmetry. 

The AdS/CFT correspondence establishes the link between the string mode solutions 
of Eqs. (2.2), (p. 3D and physical states on the gauge theory side by prescribing an UV (i.e. 



z ^ 0) boundary condition for the solutions fi (z) of the five-dimensional field equations 
Q. More specifically, for the dual of states \i) with four-dimensional spin one has to 
select the solution which behaves as fi (z) z^' where Aj is the conformal dimension 
of the lowest-dimensional gauge theory operator which creates the state The wave 
functions of states with four-dimensional spin cij acquire an extra boost factor z~'^\ so 
that the boundary condition generalizes to 

/i(z)^z^% Ti = \-ai (2.4) 

where the scaling dimension Aj of the gauge-invariant operator is replaced by its twist Tj 
Q. The lightest string modes are then associated with the leading twist operators, and 
therefore with the valence quark content of the low-spin (i.e. spin 0, 1/2, 1, and 3/2) 
hadron states [0, The duals of their orbital excitations (which have no counterparts in 



*The dependence on the four dimensions x and on the fifth dimension z factorizes at least in the 
asymptotic AdS region. 

^The Klein-Gordon, Dirac and Rarita-Schwinger equations on AdSs are discussed e.g. in Refs. [ po| , 
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the supergravity spectra) and hence of higher-spin hadrons are identified with fluctuations 
about the AdS background |l^ . 

The leading-twist interpolators contain the minimal number of quark fields q necessary 
to determine the valence Fock states of the hadrons. Their intrinsic orbital angular mo- 
mentum L is created by symmetrized (traceless) products of covariant derivatives D^. This 
results in the operators OM,f=L+2 = Q^Ds^i^ . . . Di^yq with F = 1,7^,7^ for scalar, pseu- 
doscalar and vector mesons, and Os,f=L+3 = qD^i-^ . . . Dn^qDn^^^ . . . D^^yq corresponding 
to spin-1/2 (or 3/2) baryons, with L = 'Y^^^i- The boundary condition ( p.4| ) is imposed on 
the solutions by setting the values of the five-dimensional masses m^^n (which determine 
the small- 2; behavior) according to the twist dimension f of the hadron interpolator to 
which they are dual [||, ^, 21 1, i.e. 



Y^f (r — d) for spin-0 mesons, 

\/f{T — d) + d — 1 for vector mesons, (2-5) 
f — 2 for baryons. 



where d is the dimension of the boundary spacetime. The twist dimension of the interpo- 
lating operators thereby enters the field equations and the potentials Vm,b- 

As outlined above, the duality of string modes to hadrons is based on their associa- 
tion with the lowest-dimensional, gauge-invariant QCD interpolators of matching quantum 
numbers. Hence this identification is incomplete as long as fundamental quark flavor and 
chiral symmetry are not properly accounted for. In the gravity dual, fundamental flavor 
arises from open string sectors with the strings ending on added D-brane stacks |^2| . This 
mechanism has an effective bulk description in terms of chiral gauge symmetries |9| 
whose implementation into our framework will be left to future work. 

Nevertheless, a few comments on these issues and their impact especially on the baryon 
sector may be useful already at the present stage. The simplest and currently most popular 
top-down approach to quark flavor in the meson sector introduces a D7-brane stack into 
(potentially deformed) AdS, with open strings stretching between the D7- and D3-branes 
and with the duals of flavored mesons living on their intersection |32, 33, 34]. Baryons 
require an additional D5-brane (wrapped around a compact part of the ten-dimensional 
space) which contains baryonic vertices |3^ whose attached strings pick up flavor by ending 
on the D7-branes as well [03| . 

This construction is difficult to handle explicitly, however, and the simpler propagation 
of fermionic fluctuations on the D3/D7 intersection without the D5-brane was therefore re- 
cently studied as a preparatory step |^ . The resulting modes are dual to the superpartners 
of the mesons investigated in Ref. |^3|, i.e. they correspond to fermionic bound states of 
fundamental scalars (squarks) and fundamental as well as adjoint spinors. These modes 
are solutions of Dirac equations similar to the one used in our framework but fall into 
degenerate (AA = 2) supermultiplets with their mesonic partners. 

Although bottom-up duals of baryons (with the correct 0{Nc) scaling behavior of their 
masses) obey Dirac equations as well |31, |3^, they differ from the above "mesinos" in their 
assignment to gauge-theory operators. Consequently they are subject to other AdS/CFT 
boundary conditions and generate different mass spectra, as expected on physical grounds. 



- 5 - 



Such bottom- up descriptions of baryon duals were considered in Ref. [36|, on which we 
partially rely here, and also integrated |2^] into the approach of Ref. 

Another interesting dual representation of baryons has recently emerged from a D4 /D8- 
brane construction for fundamental flavor and chiral symmetry This approach gen- 
erates a Chern-Simons gauge theory in the bulk whose instantons are dual to solitonic 



baryons of Skyrme type [37|, i.e. to collective excitations of the meson fields which carry 
topological baryon number and have masses of 0{Nc) as well. Since low-energy properties 
of Skymions can be described by Dirac fields (as they appear e.g. in chiral perturbation 



theory |3q]), this approach to baryon duals may in fact be complimentary to those of Refs. 



p5| , pi) , pq] outlined above. 

The conformal potentials induced by the pure AdSs metric^ (i.e. by Eq. (2.1) with 
A = 0) are proportional to l/z^: 



15 

— - (d - 1) 5m,v + m^uR^ 



1 



(2.6) 



(z) = m^,BR [m^,BR T 1] ^. (2.7) 



Hence the normalizable eigensolutions of Eqs. (2^) and ( |2.3D are Bessel functions whose 
order and eigenvalues depend on the boundary conditions for the solutions fi {z) of the 
field equations (cf. Ref. 0). The conformal invariance inherited from the AdS metric, 
however, prevents these potentials from carrying direct information on IR effects of QCD. 

The simplest way to approximately implement such IR effects, and in particular con- 
finement, is to impose a Dirichlet boundary condition on the string modes at a finite IR 
scale Zm- This approach is prevalent among current bottom- up models and amounts to a 
sudden onset of conformal symmetry breaking by a "hard-wall" horizon of the metric at 
the IR brane H], i.e. 

g2Ahw(.) = ^ (^^ _ ^) ^ = A^l^^^ (2.8) 

and reduces the five-dimensional, noncompact space to an AdSs slice. Even this mini- 
mal implementation of non-conformal IR effects into approximate QCD duals (with only 
one free parameter related to Aqcd) can already predict a remarkable amount of hadron 
physics, as outlined in the introduction. The investigation of hadron spectra and wave 
functions in the approach of Refs. [0, in particular, gave a good overall account of the 
angular momentum excitation spectra for both mesons and baryons. 

In view of its simplicity, however, it is not surprising that the hard wall confinement 
also reveals shortcomings. In particular, it predicts the square masses of radially and 
orbitally excited hadrons to grow quadratically with and L [^, |l^, in conflict with the 



^The modes dual to baryons originate from the ten-dimensional Dirac equation. Hence the nonvanishing 
eigenvalue of the lowest-lying Kaluza-Klein (KK) mode of the Dirac operator on the compact space X5 
| p9| adds to the five-dimensional mode mass m5,s (cf. e.g. Ref. Q). Since the AdS/CFT boundary 
conditions replace the whole mass term by a function of the twist dimension of the gauge theory operator 
to be sourced, however, the KK eigenvalue will not appear explicitly in the final expressions and has already 
been absorbed into ms b- 
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linear Regge-type trajectories found experimentally^ and expected from the semiclassical 
treatment of simple, relativistic string models [^]. While more detailed implementations 
of conformal symmetry breaking were able to resolve this problem in the meson sector 
1 11, 19, 23], linear baryon trajectories have so far not been obtained from a gravity dual. 



As mentioned in the introduction, this provides part of our motivation to search for a 
holographic representation which reproduces linear trajectories in the baryon sector as 
well. (Note, incidentally, that the approach of Ref. [11|, at least in its simplest form where 



a dilaton $ (z) oc is solely responsible for conformal symmetry breaking, will not lead 
to linear trajectories in the baryon sector since the dilaton interaction can be factored out 
of the Dirac equation and hence does not affect the baryon spectrum.) 

3. Linear trajectories of radially and orbitally excited hadrons from AdS 
type potentials 

We are now going to develop a gravity dual which manifests soft conformal symmetry 
breaking directly in the potentials and is capable of generating linear trajectories in both 
meson and baryon spectra. To this end, we find in the present section suitable potentials 
heuristically and show that they indeed reproduce the trajectories (|l . l|) . In the subsequent 
Sec. § we then construct the IR deformations of the metric (2T) which encode them 
holographically. 

A natural guess for the z dependence of potentials Wj- o which are able to generate 



the linear trajectorial (LT) structure (O) is that it should be of oscillator type in the 
infrared (i.e. quadratically rising with z for z oo). The more challenging question is 
how to realize this behavior in a universal way, i.e. on the basis of just one a priori free 
mass scale A and such that the same slope W and N + L dependence emerges in both 
meson and baryon channels. It turns out that this can be achieved at the level of the 



twist dimensions which enter the five-dimensional mass terms according to Eq. (2.5) after 
imposing the AdS/CFT boundary conditions (|2.^). Indeed, all the necessary information 
on conformal symmetry breaking can be implemented into the AdSs potentials ( p. 61 ), (^ 
by replacing 



r. + A^zl (3.1) 



(The hard wall restriction ( p^ ) of the AdS space becomes obsolete.) The heuristic rule (^ 
implies that the product of the five-dimensional masses mi and the square root a{z) = R/ z 
of the AdS warp factor grow linearly with z for z ^ oo, thus foreshadowing the linear 
trajectories ( |1 . 1] ) for both mesons and baryons. The role of the hadron-independent mass 
scale A will become more explicit below and in Sec. ^ where we relate it to the trajectory 
slope W and to the QCD scale. 

As expected from soft conformal symmetry breaking, the replacement ( |3.lD does affect 
neither the z — > behavior of the field equations nor that of their solutions. Both the 



^In the hard wall model the first radial excitations of light mesons and the nucleon, identified by a node 
in the string mode, appear at masses of about 1.8 GeV and 1.85 GeV, respectively [^], and are therefore 
difficult to reconcile with the experimental 7r(1300), p(1450) and Roper N(1440)Pii resonances. These 
shortcomings are suspected to be artifacts of the hard wall metric as well [M. 



-7- 



conformal symmetry on the UV brane and the z — > boundary conditions ( |2.4D from the 
AdS/CFT dictionary are therefore preserved. As a consequence of the above procedure, 
the mass terms in the pure AdSs potentials ( |2.6| ), ( |2.7D carry ah information not only on 
the twist dimension (and thus orbital excitation level) of the dual QCD operators but also 
on the deviations from conformal behavior in the infrared. The underlying physical picture 
will be discussed in Sec. ||. 

Recalling the expression tm = L + 2 for the twist dimension of the meson interpolators 
from Sec. |2|and making use of the replacement (3.1) then turns the mesonic AdS potential 
iK^) into 



:AV+L)^-i 



3 (3-2) 



(which holds for both spin and 1) while the AdS potential ( |2.7| ), associated with the 
baryon interpolator of twist dimension tb = L + 2>, becomes 

y]LT) (^) ^ + 1) + 1 ^ 1) + [2 (L + 1) ± 1] \^z^ + A^/} 1. (3.3) 



The normalizable solutions of the corresponding eigenvalue problems (^3) and ( |2.3| ) can 
be found analytically. For the mesons one obtains 

^^^^{z) = Nm,l,n {\z)^^"' e-^^^V2 (^2,2) (3,4) 

where the L^"* are generalized Laguerre polynomials |40| and Mh;L,n sue normalization 
constants. For the (spin 1/2 and 3/2) baryons one similarly finds 

Note that all eigenfunctions have appreciable support only over short distances, in the 
small region z < \/2A^^ ~ ^qcd i'^^- H) close to the UV brane, which is an expected 
consequence of confinement. 

The corresponding eigenvalues 

2 . . o / - - , 1 



Mil = ^X'\N + L + -y (3.7) 
Ml = A\^{n + L + (3.8) 

show that the square masses of both mesons and baryons are indeed organized into the 
observed N + L trajectories. Moreover, the spectra ( p. 7] ) and (|3.8| ) predict the universal 
slope 

W = (3.9) 

for both meson and baryon trajectories in terms of the IR scale A. They also exhibit a 
mass gap (of order VW), another hallmark of confining gauge theories, and the intercepts 
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M^Q (cf. Eq. (|1.1|)) relate the slope of the trajectories in a new way to their ground state 
masses, 



W 
3W 



The quantitative implications of these relations will be discussed in Sec. |6[ 
4. Derivation of the equivalent IR deformations of AdSs 



(3.10) 
(3.11) 



Although the existence of potentials ( |3.2| ) and (3.3) which generate linear trajectories of the 
type (|1.1| ) is encouraging, we have not yet provided any dynamical justification for them. 
Indeed, for the spectra (^7^), ( |3.8D to be the outcome of a dual gauge theory, and for the 
hadronic quantum numbers to be associated with the correct interpolating operators, one 
has to show that they emerge from stringy fluctuations in a bulk gravity background. In the 
present section we are going to establish this missing link by constructing the corresponding 
background metric explicitly. 

Of course, a priori the existence of such a bulk geometry is far from guaranteed, given 
the quasi ad-hoc nature of the heuristic rule (|3.lD which we used to find the potentials in the 
first place. Moreover, it will prove sufficient to consider just the minimal set of background 
fields^, consisting of the metric only. In fact, we will show that even the simplest type of IR 
modifications of the AdSs metric, due to a non-conformal warp factor e^^^^^^ as anticipated 
in Eq. ( |2.lD , can generate the potentials ( |3.2| ), ( |3.3| ). The success of this minimal approach 
can be at least partially understood by noting that the potentials contain effects of an order 
of magnitude which should arise from leading-order contributions to the effective gravity 
action, i.e. from the metric. Higher-order contributions due to dimensionful background 
fields (as e.g. the dilaton), in contrast, would be suppressed by potentially large mass 
scales. 

In order to prove the above assertions and to construct the non-conformal warp factor, 
we first obtain the five-dimensional field equations for string modes of the form (j) (x, z) = 
fs (z) e~^^^ (spin 0) and Vz {x, z) = fy (z) £ze~^^^ (spin 1) which are dual to mesons and 
propagate in the background of the metric (2.1) with an a priori unspecified warp function 
A {z). The ensuing bulk equations for the z dependent part of the string modes are 



9^ + 3 A' 



fsiz) = 



(4.1) 



with A' = dzA and the four-dimensional invariant square mass = -P^, as well as 

.Ax 2 



fv (z) = 0. 



(4.2) 



*An effective z dependence of the dual string mode masses may of course also arise from additional 
background fields, as for example from a Yukawa-coupled Higgs field ||41[. 
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The string modes dual to baryons can be decomposed into left- and right-handed compo- 
nents, 



i + r 



1 — 7" 



^(4) (2:), 



(4.3) 



where ^'(4)(x) satisfies the Dirac equation {i'j^df^ — M)^(4)(x) = on the four-dimensional 
(Minkowski) boundary spacetime. As a consequence, the iterated five-dimensional Dirac 
equation for reduces to 



1 



di + 4 A'--]d, + 2 A" + ^]+4 A'-- 



z . 



1 



1 



m5,BR— ] T m^,BR— [a' --\+M'^ 
z I z \ z 



1 



/±(^) =0 



(4.4) 



for the right and left handed modes with chiralities ij^f± = i/i- 

These equations can be translated into the equivalent form of Schrodinger-type eigen- 
value problems ( |2.2| ) and ( |2.3[ ) for ip (z) and V'± (z) by writing fsy (z) = [Xze~^)^^'^ (p (z) 
and f± {z) = (Aze~^)^ (z) which eliminates the first-derivative terms. The correspond- 
ing generalizations of the AdSs potentials (|2.6D , ( ^ ) are then read off from the eigenvalue 
equations as 



Vsiz) 
Vv{z) 



A" + 

2 

-A" + -A'^ 
2 



z 2 z"^ 

A' 1 1 

3— + o- 

z 2 z^ 



g2A 



o2A 



and 



^s,± (^;) = 'm^,BR- 



1 \ 

- + m^^BR— 
z I z 



(4.5) 
(4.6) 

(4.7) 



The pure AdSs potentials are contained in these expressions for ^ = 0. The AdS/CFT 
boundary condition, which relates the eigensolutions to the dual meson (baryon) operators 
of twist dimension tm = L + 2 {tb = L + 3), is imposed by adjusting the mass terms, 



ml sR^ = tm{tm - 4) = L2 _ 4^ 
ml yR^ = TAiifM - 4) + 3 = l2 - 1, 



m^nR = TB - 2 = L + 1, 



(4.8) 
(4.9) 
(4.10) 



as outlined above. Equating the general potentials (|4.5| ) - ( [4.7D to their heuristic counter- 
parts (|3.2D , ( ^ ) leaves us with differential equations for the corresponding warp functions. 
Their solutions (subject to appropriate boundary conditions), finally, determine the equiv- 



alent background metric of the form (2.1). 

As already mentioned, it is a priori uncertain whether there exists an approximate 
gravity dual whose IR deformation can reproduce a given five-dimensional potential and 
spectrum, simply because it may not result from a boundary gauge theory. In the above 
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approach this is reflected in the fact that the nonlinear, inhomogeneous differential equa- 
tions for A (z) may not have physically acceptable solutions. Our next task is therefore to 
construct and analyze the solution spaces of these differential equations for the heuristic 
potentials 



4.1 Baryon sector 

We begin our discussion in the baryon sector where the equation for A is of first order 
and hence has less and generally simpler solutions. Equating the potential ( |3.3| ), whose 
representation in terms of the background geometry we wish to construct, to the potential 
( |4.7| ) for general A results in the nonlinear, inhomogeneous differential equation 

± {zA' - 1) + le^ - [/ (/ ^ 1) + {21 ± 1) X^z^ + AS^] {le^y^ = (4.11) 

of first order (where I = L + 1 and it refers to the two baryon chiralities) whose solution 
determines the non-conformal part exp [2Ab {z)] of the equivalent warp factor. 

Remarkably, the exact (and essentially unique) solution of Eq. ( [4. 11]) subject to the 
conformal boundary condition (0) = can be found analytically and turns out to be 

Ab (z) = In (l + . (4.12) 

Note that the same solution holds for both baryon chiralities. The leading contribution to 
the non-conformal part of the warp factor at small z"^ ^ A~^ is therefore 

which has the form of the analogous warp factor exp {cz /2) used in Refs. H, ||] (together 
with a constant dilaton) to obtain a linear quark potential and a linear (mesonic) Regge 
trajectory. 

4.2 Scalar meson sector 

In the following we are going through the analogous construction for As in the spin-0 meson 



sector, which will turn out to be more multi-faceted. Equating the meson potential (3^) to 
its general-^ counterpart ( fl.5| ) produces again a nonlinear, inhomogeneous equation, but 
for As it is of second order: 

z^A" + - {zA'f - 3zA' + - (l2 - 4) (e^^ - l) - -A^z^ (x'z^ + 2L) = 0. (4.14) 
2 3 3 

In addition to As (0) = 0, its solutions require a second boundary condition which as of 
yet remains unspecified and will be determined below. This added freedom provides one of 
the reasons for the solution space in the meson sector to be larger and more diverse than 
in the baryon sector. 

In addition, the L dependence of the solutions As is more heterogeneous since the 
sign of the mass term m| ^R^ in the field equation can be either negative, zero, or positive 
(cf. Eq. ( [4.8D ). These three cases generate qualitatively different solution behaviors. 
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The positive sign corresponds to L > 2 and is associated with irrelevant gauge theory 
operators according to the renormahzation group (RG) classification. The solutions for 
L > 2 will turn out to be qualitatively similar to those in the baryon sector. The massless 
case m| gR^ = corresponds to L = 2 and to a marginal operator in the RG sense. The 
duals of the lowest orbital excitations L = 0, 1, finally, represent relevant operators and 
are tachyons^ with gR^ = —4, —3. In the following, we will discuss these three cases in 
turn. 

L = 0, 1: Due to their tachyonic nature, the L = 0,1 solutions are perhaps the most in- 
teresting ones. The negative mass term together with the specific inhomogeneity generated 
by the potential (|3.2| ) forces these solutions to develop a singularity at finite z = Zm, which 
restricts the spacetime to an AdSs slice^'^. The position and sign of these singularities 
depends on the second boundary condition for As (z). At z = this boundary condition 
may e.g. be imposed^^ on A'^' (0) for L ^ and on Ag" (0) for L = 0. The singularities 
have positive (negative) sign, i.e. A — > ±00, if A'g (0) (or Ag" (0) for L = 0) is chosen 
larger (smaller) than a critical value. 

Independently of the sign of the singularities, furthermore, the warp factor 

:^e2^«(^) >c(z„)>0 (4.15) 

{c{zm) = R'^ ex.p [2As (zm)] / z'^i^ = for negative singularities) of all solutions remains 
bounded from below for all z up to Zm where the dual spacetime ends. The above behavior 
provides a sufficient confinement criterion in the five-dimensional holographic description. 
The Wilson loop |l^] then shows an area law (since the strings are localized at Zm) and 
the gauge theory develops the expected mass gap Mmin > [^]. In fact, the hard-wall 
horizon ( |2.8D may be considered as a simple model for this type of behavior. It corresponds 
to the development of an abrupt negative singularity of A at z^. 

The origin and locus of the negative singularities can be understood quantitatively by 
obtaining a series solution for z <^ -v/2A~^ in the form 



i+E-4»w(^) 

n=l ^ ^ . 



(4.16) 



which already incorporates both boundary conditions, i.e. ^ (0) = and a second one to 
become explicit below. The coefficients An may be calculated by inserting the expansion 



^Scalar AdSs tachyons with masses satisfying the Breitenlohner-Freedman bound m?R^ > —cP/4 (which 
includes the cases we encounter here for d = 4) do not cause instabilities, as has been known for some time 

ii- 

^''in our case this is a geometric consequence of requiring the holographic dual to generate potentials 



which exhibit the linear spectral trajectories (3.7), (3.8). Alternatively, additional branes may restrict the 
fifth dimension in the IR Q . 



^^The analysis of the linearized approximation to Eq. (4.14) in App. ^ shows that As (0) = automat 



ically implies A'g (0) = 0, Ag (0) = 2LAV (^^ - 7) for L 7^ as weU as A'g (0) = A's (0) = Ag (0) = 0, 
Ag" (0) — —3 A* for L = 0, as a consequence of the inhomogeneity. 
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(I4.16D into Eq. The first three are (for L / 2) 



4L 

12^7' 



(L) 



L2 



1 



-12 
I2T5 



(4.17) 



In general, the inhomogeneity of Eq. (|4.14| ) determines the leading small-z behavior of the 
As, as revealed by the solutions ( |A.2| ) - ( [A.4| ) of the linearized equation given in App. 
For L = 0, 1, in particular, it forces the z dependence inside the logarithm to start out 
quadratically and yields 



As,L=OA {z) 



In 



^ + 1 (2A2 - Al) AS^ + O (A^^e) . (4.18) 

16 



L2-7 



(The second boundary condition can be read off from this expression.) Equation ( 4.1^ ) 
also implies that the solutions with L = 0, 1 turn negative for z > 0. At large z, on the 
other hand, the inhomogeneity rises tx and demands the modulus of As to grow as 



well. At some finite Zm the nonlinearity oc (e^"^ — l) in Eq. ( 4.14 ) will therefore become 



too negative (for As > and L < 2) to be counterbalanced by the derivative terms. To 
avoid conflict with the increasingly positive inhomogeneity, the solution then develops a 
singularity at z^- As already mentioned, the sign of the singularity depends on the second 
boundary condition. A negative singularity occurs if the slope of ^4^ near z = is smaller 
than a critical value such that the argument of the logarithm in Eq. ( 4.16| ) eventually 
reaches zero. This singularity is quantitatively reproduced by the series solution ( ^.16| ) as 
long as Zm lies inside its range of validity. 

If the slope of As at small z exceeds the critical value, on the other hand, the zero 
in the argument of the logarithm is avoided. Instead, the argument stays positive with 
increasing z, passes through a minimum and starts to increase until it reaches a positive 
pole singularity at finite z^. Singularities of this type lie outside the validity range of the 



expansion ( 4.16 ) and imply that the non-conformal part of the warp factor approaches a 
pole singularity as well. For L = 0, e.g., it takes the explicit form 



12 



X (z — ZmY 



(4.19) 



L = 2: For L = 2 the mass term in Eq. ( |4.14 ), and hence the strongest nonlinearity. 



vanishes. The initial curvature ^5/^=2(0) ~ — 4A^/3 is negative and again set by the 
inhomogeneity (cf. Eq. ( [A.3D ). Consequently, the solution can still develop a negative 
singularity if its slope close to z = remains below a critical value. For larger slopes the 
solutions turn positive and are nonsingular. Hence L = 2 corresponds to the intermediate 
case in which the solution space contains both singular solutions in which confinement 
manifests itself by compactifying the fifth dimension (similar to the L = 0,1 cases), and 
regular solutions analogous to those encountered for L > 2 (see below). Of course, the 
ensuing potential (3^) is (up to Zm) identical in both cases. 
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L > 2: For all L > 2 the main nonlinearity in Eq. ( [4.14| ) has a positive sign. Moreover, 
at z = the solutions start out with positive curvature Ag (0) = 2LA^/ (L^ — 7) (again 
dictated by the inhomogeneity, as can be seen from their linearized counterparts ( [A.4| )). 
In fact, the As^l>2 and the corresponding warp factors remain positive and nonsingular at 
all z. Typical numerical solutions for As with L = 0, 1, 2 and 3 are displayed in Fig. |l|. 
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z [X-'] 

Figure 1: Typical solutions As {z) for L = (full line, negative sign of singularity selected), L = 1 
(dotted line, positive sign of singularity selected), L = 2 (short-dashed, absence of singularity 
selected) and L — 3 (long-dashed). Note that the dual eigenmodes have significant support only 
for z < V2A"^ 



By construction, any background metric of the form ( |2.1| ), with As a solution of Eq. 
( 4.14 ), reproduces the potential ( |3.2D for all existing z. If the background space ends in 
z direction at z^, however, the potential ( p.2| ) must end there, too. This does not affect 
the (low-lying) spectrum as long as z^ ^ A^^. Hence for our purpose of maintaining the 
trajectory (|3.7D even when L = 0, 1 we choose the second boundary condition such that 
Zm becomes as large as needed^^. This selects the singularities of positive sign and implies 
that the corresponding field modes automatically satisfy Dirichlet (or Neumann) boundary 
conditions at Zm- 



/S,L=0,1 (Zm) 



XZr, 



-As{z„i) 



3/2 



fL=0,l 



0. 



(4.20) 



Even for radial excitations well beyond those currently experimentally accessible, the L = 
0, 1 part of the spectrum ( ^ ) remains therefore unaffected. (Note the Gaussian suppression 
of the eigenmodes (|3.4| ) for ^ 2A^^.) We will elaborate on this issue in Sec. |5li (Recall 
for comparison that the much slower decay of the string mode solutions in pure AdSs (Bessel 
functions) requires that boundary conditions at the hard IR wall have to be imposed by 
hand. This strongly modifies the spectrum - the masses become proportional to the zeros 
of Bessel functions Q - and generates the incorrect oc iV^, behavior which is typical 
for infinite square well or bag potentials.) 



^We have found numerical solutions of Eq. ([l.l4[) with > 6A ^ (for L — Q) 
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4.3 Vector meson sector 

The equation for the non-conformal warp factor in the spin-1 meson sector, 

_^2^" + ^^2^/2 _ 3^^/ + ^ (L2 _ 1) (e2A _ 1) _ 2^2 ^2 ^^2^2 ^ = 0, (4.21) 

is obtained by setting the meson potential (|3.2[) equal to the general-^ potential (WM). 



This equation differs from its counterpart (4.14) in the spin-0 sector only by the sign of 



the A" term and by the replacement (L^ — 4) — > (L^ — l) in the coefficient of the main 
nonlinearity. The inhomogeneity and its L dependence remain identical since they originate 
from the same bulk potential (|3.2| ) in both cases. Hence the analysis of the solution space is 



similar to that of Sec. [4.2| and will not be repeated in detail. Instead, we will just highlight 
qualitative differences in the solution behavior and discuss the numerical solutions of Eq. 
( |4.21 ) and their implications. 



The main differences between the solutions of Eqs. ( [4.14 ) and ( [4.21] ) can be rather 



directly traced to the two differences in the equations mentioned above. The modified 
coefficient of the exponential term reflects the fact that in the vector meson sector only 
L = corresponds to a tachyonic mode while the L = 1 mode is massless and all L > 1 



modes are massive (cf. Eq. (O)). This is a consequence of the additional unit of spin 
carried by the vector mesons. Hence the qualitative changes in the solution behavior 
occur around L = 1 (instead of at L = 2). The sign flip of the second-derivative term, 
furthermore, leads to changes in the sign of the solutions Ay (z) and in the development of 
singularities. Part of these modifications are determined by the solution behavior at small 
z which, as in the spin-0 case, can be obtained from the solutions to the linearized version 
of Eq. ( 4.21 ) as derived in App. 



The above qualitative expectations are corroborated by the numerical solutions (again 
subject to the conformal boundary condition Ay (0) = 0). Their perhaps most important 
novel feature is that the second boundary condition can be adapted to generate a common 
qualitative behavior for all L which is singularity- free. More specifically, every regular 
solution turns negative towards larger z where it decreases monotonically. For L > 2 the 



solutions start out with positive slope at z = (cf. Eq. ( A. 8 )) which implies that the 
decrease can set in only after passing through a shallow, positive maximum at finite z. 
While the L = solution (associated with the tachyon mode) stays regular for all choices 
of the second boundary condition, however, for L > one has the additional possibility 
of solutions which remain positive towards larger z and develop a singularity at z = 
similar to those encountered in the spin-0 sector for L = 0, 1. A set of typical solutions 
Av (z) for L = 0, 3 is plotted in Fig. ^. 

To summarize, in the vector meson channel the dual manifestation of confinement 
in terms of a dynamically compactified fifth dimension is not a necessity but rather an 
option for the higher orbital excitations (with L > 1). In contrast to the scalar sector, 
furthermore, the underlying singularities are not associated with tachyonic modes. 

We conclude this section by noting that the existence of simple gravity duals which 
reproduce the linear trajectories (|1.1|), as constructed above, provides additional support 



for the AdS/QCD program. Moreover, it establishes the basis for calculating gauge theory 
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1 2 3 4 5 

z [A-i] 



Figure 2: Typical solutions Ay (z) for L = (full line), L — 1 (dotted line), L = 2 (short-dashed, 
absence of singularity selected) and L = 3 (long-dashed, singularity selected). Recall that the dual 
eigenmodes have significant support only for z < \/2A~^. 



correlation functions and observables from the dual mode solutions ( p.4| ) - ( |3.6D on the 
basis of the AdS/CFT dictionary Q. Close to the AdS boundary, ah metrics found above 
have the same qualitative small-z behavior with a non-conformal warp factor of the form 
exp [cz"^) (except for scalar mesons with L = and vector mesons with L = 2), as deter- 
mined by the inhomogeneous terms in Eqs. ( 4.14D , ( [4.11 ) and hence directly induced by 
the underlying bulk potentials. This behavior suggests the formation of a two-dimensional 
(nonlocal) gluon condensate |45] and indicates its relevance for linear confinement (cf. Sec. 
|5.4| ). (For L = mesons the standard four-dimensional gluon condensate seems to dom- 
inate, on the other hand, which should be compared to the operator product expansion 
with renormalon-type corrections for the corresponding QCD correlators.) At large L, 
furthermore, the leading z dependence of the mesonic and baryonic warp factors becomes 
identical. Finally, an IR cutoff dual to Aqcd for the fifth dimension emerges as a necessary 
requirement for weakly orbitally (i.e. L = 0, 1) excited scalar mesons and as a possibility 
for L > 1 vector mesons. In our framework this compactification of the fifth dimension 
occurs dynamically. 



5. Discussion of the resulting holographic duals 

In the following section we elaborate on several important properties of the IR-deformed 
gravity backgrounds found above and discuss their physical significance. 

5.1 Singularities in the meson sector and linear confinement 

In several recent holographic QCD models, the asymptotically free and hence almost con- 
formal region (with at most weakly deformed AdS metric) is assumed to extend down to 
energies z~^ of the order of the QCD scale. It gets broken only in the infrared, by a rather 
abrupt onset of nonperturbative effects including condensates and confinement. 
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This scenario is indeed borne out dynamically by our results especially for the lowest 
(i.e. L = 0, 1) orbital spin-0 meson excitations. These are exactly the meson modes 
accessible in weakly curved dual supergravity backgrounds. As shown above, for them 
to lie on linear trajectories requires a singular metric of the type |17, ^] emerging in 



holographic duals of confining theories. In these cases the onset of confinement indeed 
occurs rather sudden, although not as sudden as in the extreme hard wall case (^). 
We have already traced the origin of these singularities to the tachyonic nature of the 
corresponding dual modes^^. In the vector meson sector, in contrast, singular metrics are 
possible but not mandatory (depending on the choice of a boundary condition) for orbitally 
excited resonances, and they are not tachyon-induced. 

Physically, it seems that no (L = 0) or only a small (L = 1) centrifugal barrier allows 
the corresponding spin-0 meson states to probe different aspects of the IR region and hence 
to feel a more sudden onset of confinement, in particular at large N. The location of the 
singularities, however, is set by the inverse mass scale in the second boundary condition and 
can therefore be put at ^ large enough not to significantly alter the low- lying part 
of the spectra. Furthermore, semiclassical arguments indicate that highly excited hadrons 
generally become larger and therefore should be able to explore more of the IR region as 
well. This may explain why the vector meson metric can become singular for L > 1. 

The possibility to set Zm J?> implies, in particular, that one could choose extensions 
of the metric into the region z G [zm — oo] which yield the same low-lying spectra and 
wave functions (cf. Eq. ( [4.20D ) without any singularities. In the intermediate case of spin-0 
mesons with L = the choice between singular and regular metrics exists even inside the 
solution space of Eq. ( [4.14 ). For scalar meson excitations with L > 2, vector mesons under 



suitable boundary conditions and all baryons, finally, confinement does not manifest itself 
in metric singularities at all. 

Nevertheless, all higher meson and all baryon excitations are found to lie on the linear 



trajectories^^ (3^) ^-nd ( |3.8D . The combination of these results may reflect the fact that the 
wave functions of light hadrons seem to be rather weakly affected by the linear confinement 
force. Indeed, several successful models for low-lying hadrons (e.g. models of Skyrme-type 
1 37] and the instanton-based chiral quark model do not implement confinement at 
all since the mean separation among colored constituents appears to be too small for 
confinement effects to become relevant. 

As noted in Ref. 0, the abrupt hard wall singularity of the metric ( |2.8D - and the 
IR boundary conditions it requires - resemble those of an MIT bag model with a sharp 
surface [|5^. In our case the singularities develop more gradually. The analogous bags 
therefore have smooth transition regions as they emerge dynamically in soliton bag models 



of Lee-Friedberg |51] or color-dielectric |52| types. Such soliton bags become confining by 



means of a space-dependent color dielectric function which induces singular couplings to 



13 



It may be tempting to speculate about potential relations between these tachyon-induced singularities 



and closed string tachyon condensation 
gauge theory side (see also Refs. 4J 



4^ in the bulk, which is expected to be dual to confinement on the 



The holographic models of Refs. [1^, |l^] also realize linear meson trajectories with a non-singular 
metric, but they contain an additional (regular) dilaton field. 
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vacuum fields |^^. This suggests that the IR deformations of the dual gravity background 
found above encode information on the color-dielectric QCD vacuum structure. 

String breaking due to light quark production is expected to stop the linear rise of 
the QCD confinement potential at large distances and hence to bend the linear hadron 
trajectories at sufficiently high excitation levels. Although such effects are not yet visible 
in the experimentally accessible part of the hadron spectrum, string breaking has recently 
been confirmed on the lattice [^]. In our model, however, the linearity of the spectra 
( |3.7| ), ( p^ ) continues up to arbitrarily high excitation quanta (except for spin-0 mesons in 
L = 0, 1 states with — > oo, due to the finite z effects discussed above). This indicates 
that string breaking effects are absent in our holographic dual, as expected in the large- 
limit (where is the number of colors) or the associated weak coupling approximation on 
the gravity side. 

A simple way to account for string breaking effects by hand would be to modify the 
solutions A at z > Zm — £ such that the potential levels off. As already alluded to, for 
(Xzm)'^ ^ 1 any reasonably smooth deformation of A at large z may in fact be implemented 
with practically no impact on the low-lying part of wavefunctions and spectra. In this way 
one could for example remove the singularities altogether. Not surprisingly, this also implies 
that the generation of the linear trajectories ( |1.1| ) does not fully constrain the IR behavior 
of the gravity background. Together with the second boundary condition in the meson 
sector, the remaining freedom could be used to implement a more comprehensive set of 
QCD observables. 



5.2 L dependence 



The IR deformations of the AdSs metric obtained from the solutions of Eqs. (4.11), ( 4.14| ) 



and ( [4.21| ) are necessarily L dependent. This L dependence enters through the potentials 
( |3.2| ), (|3.3D which give rise to the inhomogeneities of the differential equations for Asy^B, 
and more universally through their counterparts (4^) - ( [4. 7]) in the IR deformed back- 



ground. Its ultimate source is therefore the L-dependent twist dimension of the considered 
hadron interpolators, imposed via Eqs. ( ^4.8]) - ( [4.10 ), and its main effects are indepen- 



dent of the specific choice for the heuristic potentials or the replacement rule (3.1). A 
somewhat analogous hadron dependence of a dual background has been found in Ref. pO| , 
where vector and axial vector mesons feel a different metric, and would enter several other 
holographic models if observables in the whole hadron spectrum were to be reproduced. 

A natural source for the L dependence arises in our approach from the identification 
of orbital hadron excitations as stringy quantum fluctuations about the AdS background 
1^, 16 1 . Indeed, such L dependent fluctuations may deform the AdS background met- 



ric in an L dependent fashion. In the simpler case of two-dimensional quantum gravity, 
fluctuation- induced deformations of a background metric (essentially AdS2) were recently 
found explicitly [^]. In our case, the back-reaction of the metric to a fluctuation dual to a 
given orbital excitation could conceivably lead to analogous, L dependent deformations, as 
found in our solutions Asy^B (z)- Although the different orbital excitations feel a different 
total metric, however, the overall conformal symmetry breaking scale A remains (almost) 
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hadron independent. This is a consequence of its relation ( |3.9D to the almost universal 
slope of all hadron trajectories, which we will discuss quantitatively in Sec. ^. 

The identification of orbital excitations as duals of metric fluctuations sets their holo- 
graphic origin apart from that of the radial (i.e. A^) excitations. This becomes manifest 
in the fact that our metric acquires no N dependence while even more sophisticated IR 
deformations would necessarily be L dependent since each orbitally excited hadron state 
is created by a different interpolator. Via back-reactions similar to those considered in 



Refs. [48, 57] this L dependence may carry over to additional background fields of stringy 
origin as well. The relative weakness of the warp factors' L dependence in the physically 
dominant region z < \/2A^^ and for larger L (where it becomes identical in the meson and 
baryon sectors) may be another indication for their fluctuation-induced origin. 

Additional support for the above interpretation of the angular-momentum dependence 
in our background arises from an observation in Ref. The latter asserts that for on- 

shell gauge theory properties which are described by the quadratic part of the dual string 
action, like the trajectories ( |1.1| ), the effect of higher derivative terms (including those 
related to orbital angular momentum operators, cf. Sec. ^) can be essentially reproduced 
by the standard quadratic terms - to which we restrict ourselves here - in a modified 
gravity background. The modifications will depend, in particular, on the hadronic angular 
momentum carried by the dual modes, as manifested in our case in the solutions Asy^B- 

Pursuing the above line of reasoning farther, we recall that in the holographic model 
based on the hard-wall metric (|2.8| ) quantum fluctuations corresponding to orbital exci- 
tations are (at least in the conformal regime at small z) represented by an L dependent 
eff^ective mass for the bulk string modes. Our derivation of the dual gravity backgrounds 
suggests a generalization of this interpretation. By allowing the effective masses to change 
with resolution r = R'^/z outside of the conformal regime (i.e. for 2; > 0, as via the re- 
placement (p.ip) one may be describing the quantum fluctuations of the metric and their 
potentially deforming back-reaction in more detail, and hence obtain a more accurate de- 



scription of IR properties (including the linear hadron trajectories (1.1)) on the gauge 
theory side. 

5.3 Background field content and underlying dynamics 

The construction of our background geometry raises the question how it may be related to 
an underlying string theory. One could start to gain insight into this matter by examining, 
for example, whether the differential equations ( [4.11| ), ( 4.14 ) and ( [4.21[) for the metric can 



be at least approximately cast into the form of (potentially higher-dimensional) Einstein 
equations. Such issues are beyond the scope of the present paper, however, where we focus 
on general principles, symmetries and experimental data to constrain the holographic dual 
in bottom-up fashion but leave the underlying dynamics at least a priori undetermined. 
In the remainder of this section we will therefore only mention a few qualitative aspects 
of the dynamics expected to govern the dual background and summarize the rationale for 
restricting our construction to the metric. 

As stated in the introduction, the present knowledge of string theory in strongly curved 
spacetimes does not provide ab initio insight into the background fleld composition and 
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dynamics of the QCD dual. (Guidance from the supergravity approximation, in particu- 
lar, is limited even at large Nc- it would predict a mass gap far smaller than the string 
tension beyond which linear trajectories can emerge, for example, in conflict with QCD 
phenomenology p!l|].) Nevertheless, several approximate, QCD-like duals, partly including 
fundamental flavor as mentioned in Sees. |l] and ^, have recently been obtained from brane 
constructions in which the metric and other background fields are solutions of classical field 
equations. Other dual models were derived by solving simple, string-inspired scalar and 



Einstein equations |^, 57 1 and used to get an idea of the impact of QCD condensates on 
the background. These models show, in particular, that such back-reactions can generate 
confining IR cutoffs in the fifth dimension similar to those encountered in our approach. 

Our exclusive reliance on the dual background metric, finally, was guided by the di- 
mensional power-counting argument of Sec. ^ the possibility to represent the impact of 
additional background fields on the resonance masses by deformations of the metric (cf. 



Sec. 5.2), and by an Occam-razor type preference for the minimal and hence most efficient 
background required to reach our objectives. Nevertheless, experience from brane models 
as discussed above suggests that a more comprehensive and detailed approximation to the 
QCD dual should contain additional background fields. As pointed out in Ref. |11], for 
example, one may expect a background including tachyon and dilaton fields if the dual 
confinement mechanism has its origin in closed-string tachyon condensation (see also Ref. 

Hi)- 



5.4 Comparison with other confining holographic models 

In the following we will briefiy compare our holographic model to a few related approaches 
which also contain dual representations of linear confinement and linear trajectories in the 
hadron spectrum. As already mentioned, it turns out to be a nontrivial task to reproduce 
linear trajectories with approximately universal slopes not only in the meson but also in 
the baryon channels. In fact, our approach seems to be the first which accomplishes this. 
Hence our comparisons below have to remain restricted to the meson sector. 

Linear "meson" trajectories in more or less QCD-like gauge theories were e.g. found 



in Refs. |23, 42]. A recent implementation of linear trajectories for both radial and spin 



excitations of the rho meson into the AdS/QCD framework |11| induces conformal sym- 
metry breaking mainly by a dilaton background field. In the simplest case a dilaton of 
the form <I> (z) oc is added to the pure AdS metric and hence is exclusively responsible 
for the non-conformal IR behavior. (Essentially the same term was argued to arise from a 
dual magnetic condensate (which plays the role of a Higgs field) in the partition function 
of a QCD instanton ensemble when promoted into the bulk |5^.) Although this approach 
works well in the (vector) meson sector, we have already noted that such dilaton effects 
do not manifest themselves in the baryon spectrum since they can be absorbed into the 
eigenmodes and leave the AdS/CFT boundary condition unchanged. An interesting obser- 
vation of Ref. [11 1 is that the exponent of the warp factor should not contain contributions 
growing as z'^ for z — > oo, in order to have spin-independent slopes of the radial rho meson 
excitation trajectories. Our solutions Asy^B grow logarithmically with z for large^^ z and 
except if the large-z region is cut off by a singularity, of course 
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are therefore consistent with this condition (while the model of Refs. [19, 42] is not, see 
below) . 

In other recent work, scalar bulk fields dual to the gluon and bilinear quark condensate 
operators have been shown to generate - by their back-reaction on the gravity background 
- a confining restriction of the metric to a deformed AdS slice |4^, |5^ . Indications for the 
potential role of the two-dimensional QCD condensate [^] in the confinement mechanism 
arise in these frameworks as well. As already mentioned, the non-conformal behavior of our 
metric towards small z, induced by the confinement-generating modifications of the string 
mode potentials, may similarly be a reflection of the non-local dimension-two condensate. 
This suggest a more general analysis of the information on QCD condensates and more 
specific IR degrees of freedom (e.g. topological ones related to instantons, monopoles or 
center vortices) |58] which is encoded in our background metric. 

The qualitative small-z behavior exp (cz^) of our solutions for the non-conformal warp 
factors is identical to that proposed in Refs. ||l^, 42| where it has been shown to embody, 
together with a constant dilaton, a linearly growing heavy-quark potential and a linear 
(mesonic) Regge trajectory |4^. The value of c was estimated to be c ~ —0.9 GeV^ (in 
a metric with Lorentzian signature) which is in the same ballpark as ours for intermediate 
L. Indeed, anticipating the relation ( |6.5D between A and Aqcd to be established in Sec. 
|6|, our warp factor implies e.g. cs,l=2 — —0.7 GeV^. Our values for cb in the baryon 
sector are positive, however, which may suggest some differences in the dual confinement 
mechanisms for vector mesons and baryons. 



6. Phenomenological implications 



In the following section we proceed to the quantitative analysis of our holographic dual and 
confront the predicted mass spectra ( |3.7| ), ( p.8| ) for the light hadrons with experimental 
data. Recent reviews of excited hadrons, their symmetry structure, parity doubling etc. 
can be found in Refs. |59|. 

We start by determining the conformal symmetry breaking scale A from data for the 
slope W = 4A^ of the trajectories ( |1.1D . Fits to the experimental meson spectra yield 
W = (1.25 ±0.15) GeV^ ^andW = (1.14 ±0.013) GeV^ [g§. These values allow for an 
immediate check of our relation ( 3.101) which predicts the rho meson mass as a function of 
its trajectory slope, i.e. 

fw 

(6.1) 



p — 0.76 GeV, respectively, 
0.7755 ± 0.0004 GeV M- 



The above empirical results for W imply Mp = 0.79 GeV or 
which are both consistent with the experimental value Mp 
Since the latter is close to the mean of the slope fit results, we choose the experimental rho 
mass to set the scale of the deformed gravity background, i.e. 



'E = ^4 = 0.55 GeV. 

4 ^/2 



(6.2) 
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The corresponding value W = 4A^ = 1.21 GeV^ fixes the slope of our linear meson trajec- 
tory which is compared to the experimental meson resonance spectrum (for quark-antiquark 
states) in Fig. ^. 
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Figure 3: Experimental meson mass spectrum from Ref. 
W = 2M2 ~ 1.01 GeV^. 



pOj] and the predicted trajectory for 



The clustering of radial and orbital excitations is clearly visible, and even the highest 
radial excitations /2(2300) and /2(2340) lie squarely on the linear trajectory. The pion 
ground state, set apart by its approximate Goldstone boson nature, does not fit into the 
overall pattern predicted by the dual string modes. This problem is expected. It was 
already encountered in Ref. ||7| and is caused by the lack of chiral symmetry and its 
breaking in our approximate holographic dual. For the same reasons, qualitative models 
for the light-front square mass operator in the mesonic valence quark sector (which 
reproduce the radial excitation trajectory) need a strong additional short-range attraction 
in the spin-0 channel to reproduce the n-p mass splitting (for L = 0). Since the dual 
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string modes are related to the valence components of the light-front wave function (with z 
playing the role of a relative coordinate) ||2l| , the impact of such interactions is accessible 
in our approach. 

The empirical slope of the A trajectory (including the nucleon resonances in the ^8 
representation of SU(4)) is = (1.081 it 0.035) GeV^ ||2^. As in the meson sector, our 
relation ( p. 11 ) turns this value into a prediction for the A ground state mass, 



Ma = -J — = 1.27 GeV, (6.3) 

which compares well with the experimental value Ma = 1.232 GeV. This suggests to use 
the experimental mass of the A isobar for an alternative determination of the scale 

A = ^ = 0.50 GeV (6.4) 
V6 



in the baryon sector. The value ( |6.4| ) differs by less than 10% from that based on the 
experimental rho mass, Eq. (|6.2|). This confirms the approximate universality of A and of 



the associated slopes W = 4A^ in the meson and baryon sectors. 

The resulting A isobar trajectory, together with the empirical square masses of the 
radial and orbital resonances, is shown in Fig. ^. The first radial excitations, A(1600)|^ 
(with L = 0) and A(1930)| (with L = 1) are degenerate with states carrying one or two 
units of angular momentum, respectively. The parity doublets A(1600)|^, A(1700)| and 
A(1905)|''', A(1930)| , furthermore, are states which differ by one radial and one orbital 
excitation quantum such that + L is preserved. 

Finally, we turn to the nucleon and its excitations. As noted in Ref. the trajectory 
of the nucleon resonances in the ^8 representation of SU(4) (including the nucleon itself) 
lies below that for the ^8 representation. This behavior can be accommodated by our 
covariant^^ framework with a somewhat smaller value of A = 0.47 GeV which may e.g. be 
due to hyperfine interactions. The resulting trajectory indeed fits the nucleon resonances 
in the ^8 representation well and is shown as a solid line in Fig. |5[ The experimental value 
of the nucleon mass lies below the value Mjv = 1.16 GeV on the trajectory, however. This 
may again be related to chiral symmetry breaking effects and implies, in any case, that the 
resulting nucleon delta splitting vanishes. Since the latter is an O (l/Nc) effect (where Nc is 
the number of colors), this result is consistent with the large- Ac limit, i.e. the weak string 
coupling limit which underlies all known top-down holographic duals. The trajectory of 
Fig. Q is included as a dashed line in Fig. |5| and seen to fit the nucleon resonances in the 
^8 representation, as anticipated. 

Our holographic results for the at present experimentally accessible orbitally excited 
nucleon states are generally close to those of Ref. [^] (which were based on the hard IR 
wall metric ( |2.8| )) because at moderate N + L the difference between linear and quadratic 
trajectories is rather small. In addition, we predict the radially excited states, i.e. the Roper 
resonance A(1440) and the second radially excitation A(1710), which is almost degenerate 



^^Note that the identification of the not separately Lorentz-invariant orbital angular momentum L requires 
us to select a particular frame, which underlies the interpretation of our interpolators (cf. Sec. H). 



-23- 



8 - 



6 - 



> 



4 - 



2 - 







A(2420) 



A(1700) 
A(1620) 
A*(1600) 



A(1950) 

A*(1930) 
A(1920) 

A(1910) 

Al905) 




A(1232) 







2 3 
N+L 



6 



Figure 4: Experimental Delta isobar mass spectrum from Ref. | 
for S = 3/2 (in the ^8 representation of SU(4)) with W = 2MI/3 



I and the predicted trajectory 
1.01 GeV^. 



with the L = 2 states A^(1680) and A^(1720). Our results for nucleons with internal spin 
3/2 lie approximately on the A trajectory (dashed line), as already mentioned. The parity 
doubling of baryon states with fixed total spin, differing by one unit of angular momentum 
and one internal spin or radial excitation quantum, emerges naturally in our approach. 

The slope W of our trajectories is related to the QCD scale. As a consequence of 
confinement, the Gaussian suppression factor exp [— {W/8) z^] prevents the dual string 



modes (3^) - ( |3.6D from extending significantly beyond distances Zm ~ y/sJW into the 
fifth dimension. Similar confinement effects are often modelled by the hard IR wall metric 
( |2.8| ) with Zm ~ ^QCD' ^ discussed in Sec. Hence we approximately identify 

W A 



Aonn ~ A / — = ~ 0.35 GeV. (6.5) 
The numerical estimate in Eq. ( |6.5| ) is based on the phenomenological slopes of about 
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Figure 5: Experimental nuclcon mass spectrum from Ref. ||60| and the predicted trajectories for 
S = 1/2 with W ~ 0.9 GeV^ (solid line) and for S = 3/2 with W = 2APji ~ 1.01 GeV^ (dashed 
line). 



1 GeV^ and indeed close to the empirical value Aqcd — 0.33 GeV (at hadronic scales 
with three active flavors) |^0|. Finally, we recall that the string tension resulting from the 
semiclassical treatment of simple, relativistically rotating string models |2^, i.e. 



ivr fm 



is consistent with standard values as well. 



7. Summary and conclusions 

We have shown how a salient empirical pattern in the light hadron spectrum, namely the 
combination of both radial and orbital excitations into linear trajectories of approximately 
universal slope, can be reproduced with good accuracy by a rather minimal version of 
holographic QCD. Our approximate holographic dual relies exclusively on IR deformations 
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of the AdS metric, governed by one free mass scale proportional to Aqcd, and generates the 
mass gap expected from confining gauge theories. Moreover, it provides the first example 
of a gravity dual which is able to reproduce linear trajectories in the baryon sector as well. 

The resulting light hadron spectra are in good overall agreement with the available 
experimental data for both meson and baryon masses. Discrepancies between the radial 
and orbital resonance masses in the hard-wall model are resolved and the experimentally 
established, approximately universal slope of the trajectories emerges naturally. Moreover, 
new relations between the p meson and A isobar ground state masses and the slopes of their 
respective trajectories are predicted. Since the linearity of all trajectories extends to the 
lightest masses, however, they fail to reproduce the physical pion and nucleon ground states. 
This is not unexpected because our approximate gravity dual in its present form lacks 
information on chiral symmetry and residual interactions responsible e.g. for hyperfine 
splittings. 

Our holographic background was derived by reconstructing the dual mode dynamics 
from spectral properties on the gauge theory side. The underlying strategy may be useful 
for other applications as well. It consists of first finding the modifications of the AdS 
string mode potentials which generate a desired gauge theory result, and to subsequently 
construct the corresponding fields on the gravity side by equating the potentials induced 
by a general background to their heuristic counterparts. As long as a dual background 
exists, its derivation is then reduced to solving the resulting differential equations. 

Above we have constructed the holographic duals for radial and orbital hadron trajec- 
tories in a minimal way, i.e. by a non-conformal warp factor. A remarkable a posteriori 
justification for this restriction is its sufficiency. More complex IR deformations and fur- 
ther bulk fields would also introduce new parameters to be fixed by QCD phenomenology 
and hence lessen the prcdictivity of the dual description. In the baryon sector we were 
able to derive the resulting IR deformations of the metric analytically. The lowest orbital 
excitations of the spin-0 mesons, and depending on a boundary condition also the higher 
orbital vector meson excitations, encounter a singular metric. This is a dual confinement 
signature and results in a dynamical compactification of the fifth dimension, hence directly 
linking linear trajectories (for not too high excitation levels) to linear quark confinement. 

Despite the advantages of the minimal description, however, experience from super- 
gravity and brane models suggests that a more comprehensive holographic dual may require 
a more general form of the metric and additional background fields. The prospect of de- 
riving those by our method deserves further investigation. Among the potentially useful 
extensions and applications we mention the implementation of quark flavor and sponta- 
neously broken chiral symmetry as well as the calculation of condensates, heavy-quark 
potentials and light-front wave functions. 
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A. Solutions of the linearized differential equations for the mesonic warp 
factor 



In this appendix we solve the linearized versions of the differential equations (4.14) and 



(4.21) for the non-conformal mesonic warp factors. The results will be useful for our analysis 



of the solution behavior of the full, nonlinear equations in Sec. ^ It will also be instructive 
to understand the differences in the solution behavior induced by tachyonic, massless and 
massive modes in this simplified setting, although the solutions of the linearized equations 
do not develop finite-z singularities. 

The linearization of the differential equation ( [4.14 ) for the spin-0 meson sector leads 
to the (still inhomogeneous) equation 

z^A" - 3zA' + 1 - 4) A - h'^z^ (A^z^ + 2L) = (A.l) 

whose solutions provide approximations to those of the full equation in the regions where 
A <^ 1. As a consequence of the conformal boundary condition ^(0) = 0, which turns 
out to be an automatic property of all solutions which stay finite at z = 0, this condition 
should hold in particular in the UV, i.e. for z close to zero. 



The full solution space of the linear equation (AJ) can be constructed by Frobenius 
expansion techniques (for the homogeneous part) and by guessing special solutions of the 
full, inhomogeneous equation or by deriving them with the help of Green function methods. 
Either way, the general solution for L = 0, 1 (associated with the tachyonic string modes 
in the bulk) is found to be 

v2 9 x4^4 



As,L=OA (z) = + ^^3^ + ci (Az)2+V— + C2 (Az)2-V— . (A.2) 

This solution satisfies the initial condition ^ (0) = only for C2 = 0. The remaining 
irrational-power term is subleading at small z both for L = and L = 1. For L = 2, the 
general solution (induced by the massless string mode) is 

As,L=2 {z) = -h'z^ - 7^ (1 + ci) X^z" + Ix'zHnXz + C2 (A.3) 
3 24 

where the initial condition A{0) =0 requires C2 = 0. For L > 2, finally, the particular 
solution of the inhomogeneous equation is identical to that for L = 0,1 while the general 
solution of the homogeneous equation differs. Their sum, 

LX^z^ AS4 



As,L>2 (z) = -p^ — t; + 



L2-7 2L2 



+ ci A^z^ cos 2\ — - — InXz] + cgA^z^ sin 2\ — - — In Az , (A.4) 
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is the general solution of equation ( A.l ) and satisfies the initial condition A{0) = for 
all (finite) values of ci^2- Note that for L = 0,1 the coefficients of both X^z^ and X'^z^ 
are negative, for L = 2 the coefficient of A^z^ is negative and that of A^z^ can be chosen 
positive (so that Ai=q > for larger z), and for L > 2 both coefficients are positive. 
The analysis of the linearized version 



-z^A" - 3zA' + 1 (l2 - 1) ^ - -X^z^ {X^z^ + 2L)=0 



(A.5) 



of equation (4.21) in the vector meson channel proceeds analogously. For L = it yields 
the general solution 



Av,L=o {z) = -\\'z' - ^AS^ + ci (Xz)-' cos 



, . /'InAz 

C2 (Az) sm ( 



where Av,l=o (0) = demands ci = C2 = 0. For L = 1 one finds 

Av,L=l (z) = - Ja^Z^ _ Lx^z^ + ci-1- + C2 

3 36 X z^ 

where Av,l=i (0) = again requires ci = C2 = 0. The solution for L > 1, finally, is 



W,L>1 



{z) 



LX'z 



2^2 



4^4 



+ 



A*z 



L2-7 2L2-38 



+ ci (Az)" 



+ C2 (Az) 



(A.6) 



(A.7) 



(A.^ 



where Av,l>i (0) = demands ci = while C2 remains unconstrained. For L = 2 the 
irrational power term provides the leading small-z behavior. The coefficient of the A^z^ 
(A'^z^) term turns positive for L > 3 (L > 4). 

determine the small-z behavior of the so- 



The inhomogeneities in Eqs. (A^) and (^ 
lutions (with the exception of Av,l=2)- In our context, this has two pertinent consequences. 
First, the leading small-z dependence is generally determined by the special solutions of the 
inhomogeneous equation and therefore completely fixed, i.e. only the subleading small-z 
behavior depends on the boundary conditions. 

A second useful consequence of the inhomogeneities in Eqs. ( [A.l| ) and ( [A.5| ) - which 



are the same as those in Eqs. (4.14) and (4.21) - is related to the fact that the leading 
small-z; behavior of the solutions to the full equations is identical to that of their linearized 
counterparts. Hence at small z the modulus of all solutions with ^ (0) =0 grows as A^z^ 
(except for As^l=o and Av,l=2), partially with trigonometric or logarithmic corrections. 
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